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The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma 
Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRBs) 
is isotropic but is bound radially (Meegan et a!. 1992). This finding suggests that the bursts 
are either cosmological (Hartmann 1992; Paczynski 1994) or they originate from an extended 
Galactic halo (Eichler & Silk 1992; Li & Dermer 1992). The implied luminosities and the 
observed variability of the GRBs on time scales as short as one millisecond suggest that they 
originate from compact objects. 

We are presently studying black hole models for GRBs. Any such model must produce 
a non-thermal photon spectrum to agree with the observed properties. For a wide range of 
burst parameters the assumed bursting source consists of a non-thermal electron-positron- 
photon plasma of very high density. It seems possible to produce such a plasma in accretion 
onto black holes (Gunnlaugur & Svensson 1992). 

Early attempts to understand the physical properties of pair plasmas were restricted 
to the case where the pairs were assumed to be in thermodynamic equilibrium. Bisnovatyi- 
Kogan, Zel’dovich, & Sunyaev (1971) analyzed pair plasmas, with negligible photon concen- 
tration, at extreme relativistic temperatures. Their work was generalized and enlarged by 
Lightman & Band (1981), Lightman (1982), and Svensson (1982), who considered thermal 
pair plasmas in equilibrium with photons. These latter papers contain results for a wide 
range of temperatures and include many novel photon and particle reactions, which were 
comprehensively tackled for the first time. They identify various parameters of a thermal 
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pair-photon plasma and give a detailed portrait of plasma equilibria in the parameter space. 
All these papers assume that the plasma is stationary and homogeneous. Guilbert & Step- 
ney (1985) and Kusunose (1987) analyzed the time evolution and radiative-transfer problems 
in a thermal pair plasma with the aid of kinetic equations. Investigation of the physics of 
non-equilibrium pair plasmas and their time evolution is a relatively recent undertaking. 
Nonetheless, several important results have been obtained already. Those studies, in gen- 
eral, have employed two different methods: Monte-Carlo method (e.g., Stern et al. 1995 and 
references therein) and solving the kinetic equations (e.g., Coppi 1992; Lightman & Zdziarski 
1987; Svensson 1987; Zdziarski, Coppi, & Lamb 1990). 

In our on-going work, we are developing the kinetic theory for a non-equilibrium pair 
plasma. For such a plasma, containing arbitrary concentrations of pairs and energetic pho- 
tons, we derive the expressions for the emission and absorption coefficients, for particles 
as well as photons, starting from a single general equation. The processes we consider are 
Compton scattering, pair annihilation and creation, and bremsstrahlung. The kinetic kernels 
for all these processes are given in terms of arbitrary distribution functions, for pairs as well 
as photons. We then use this formalism to compute the luminosities, cooling rates/opacities, 
and the relaxation times for such a plasma. 

The main new features of our work are as follows. (1) We do not assume the presence 
of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an 
arbitrary concentration and energy distribution. This means that the kernels for processes 
such as Coulomb collisions are not linear functionals of the distribution function any more. 
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This situation affects the photon distribution also, in a non-trivial way, because the kinetic 
equations are coupled and are non-linear. (3) There is no soft photon source in our model. 
Initially all the photons in the plasma are either energetic X-rays or 7 -rays. Therefore most 
of the photons are initially very near or above the pair-creation threshold energy. Also, there 
are no cold electrons in the plasma. (4) The initial energy distribution of the pairs as well as 
photons is arbitrary. Also, the pair and photon concentrations may be comparable to each 
other in the beginning. (5) We collect the analytical expressions for the kinetic kernels for 
all relevant processes. Several of these are evaluated through a method that is different from 
those that are available now. Where possible, we compare our results with previous ones. 
( 6 ) We present a different approach to finding the time-evolution of pair and photon spectra, 
which is a combination of the kinetic-theory and the non-linear-Monte-Carlo schemes. 

We have developed many Monte-Carlo programs to model various processes, to take 
into account the time evolution, and to incorporate various physical effects which are unique 
to non-thermal plasmas. We have tested our programs by compaxing our results with those 
of the existing ones for the special case of equilibrium distributions and they are in very 
good agreement. All results pertaining to this study, in full detail, are in a paper (Pilla and 
Shaham 1995) in the final stages of preparation. 

Fireballs are a generic feature of many gamma-ray-burst models. The hydrodynamics 
of fireballs in GRBs was studied before (e.g., Shemi & Piran 1990; Piran, Shemi, & Narayan 
1993; Meszaros, Laguna, & Rees 1993; Meszaros & Rees 1993). Applying results from kinetic 
theory will improve our understanding of these systems. While the full understanding of the 
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nature of non-thermal radiation from various compact sources, especially GRBs, requires 
the solution to the radiative transfer equations which incorporate Comptonization, pair 
processes, and bremsstrahlung in a self-consistent manner, the analysis carried out so far 
serves as a preliminary to such studies. In a forthcoming paper (Pilla 1995) we will apply 
the formalism of the present study to the fireballs in GRBs. 

The research described above has been carried out by R. Pilla with the late Jacob Sha- 
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